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Frozen states and active-absorbing phase transitions of the Ising model on networks 
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A zero temperature quench of the Ising model is known to lead to a frozen steady state on 
random and small world networks. We study such quenches on random scale free networks (RSF) 
and compare the scenario with that in the Barabasi-Albert network (BA) and the Watts Strogatz 
(WS) addition type network. While frozen states are present in all the cases, the RSF shows an 
order-disorder phase transition of mean field nature as in the WS model as well as the existence of 
two absorbing phases separated by an active phase. The WS network also shows an active-absorbing 
(A-A) phase transition occurring at the known order-disorder transition point. The comparison of 
the RSF and the BA network results show interesting difference in finite size dependence. 
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I. INTRODUCTION 

Dynamics on networks is a topic on which extensive 
works have been done in recent years. Phenomena which 
have been studied include evolution of spin systems, opin¬ 
ion dynamics, disease spreading dynamics, etc. 0 , 0 . 
The dynamical picture is quite different from that on reg¬ 
ular lattices due to the topological features of network. 
For example, one can define the dynamics in different 
ways for the voter model on a network while these rules 
become identical on lattices 0. 

The study of Ising model on networks has revealed a 
number of interesting features when static properties are 
considered. Even in one dimension, when randomly new 
links are added (or existing links rewired) as in a Watts 
Strogatz (WS) network 0, one gets a phase transition 
which occurs with mean field criticality [HI3. On 
Euclidean networks, indications of both mean field type 
and finite dimensional-like phase transitions have been 
shown to exist by var ying th e relevant parameter 0. On 
scale free networks [12 h 14| . the transition temperature 
shows a logarithmic increase with the system size which 
is perhaps the most surprising result [15h4l8l | . 

While considering ordering dynamics on regular lat¬ 
tices at zero temperature for the Ising model using 
Glauber dynamics, it is known that for any dimension 
greater than one, freezing occurs with a probability de¬ 
pendent on the dimension 0 . This happens when one 
considers a completely random initial condition. On ran¬ 
dom graphs and networks, one encounters similar freezing 
phenomena which depend on the density of added links 

On random networks or graphs, the evolution of the 
Ising model from a completely random state shows that 
the system does not order. A freezing effect was observed 
and although there could be an emergent majority of 
nodes with either spin up or spin down state, domains 
of nodes with opposing spins survive HKHl. Careful 
observations show that the disordered state is not an ab¬ 
sorbing state 0 • It is instead a stationary active state, 
with some spins flipping, while keeping the energy con¬ 
stant. The number of domains remaining in the system 


is just two. The qualitative picture is then the same as 
on regular lattices for d > 2 [l{|, the system wanders 
forever in an iso-energy set of states. The distribution 
of the steady state residual energy (which is identical to 
the number of bonds between oppositely oriented spins 
apart from a constant) for the Ising model on a random 
networks was investigated j2{|. It was found that the 
distribution typically shows two peaks, one very close to 
the actual ground state where the residual energy is zero 
and one far away from it. 

Dynamics of the Ising model on the WS network with 
restricted rewiring has also been considered. Here ini¬ 
tially a spin is connected to its four nearest neighbours 
and then only the second nearest neighbour links are 
rewired with probability p. The system therefore always 
remains connected. Under the zero temperature Glauber 
dynamics, freezing effect was observed for any p ^ 0 [26| . 

In this paper, we have considered in detail the variation 
of the relevant thermodynamic quantities as functions of 
time for the zero temperature dynamics of Ising model 
on various networks. Our main emphasis is on the Ising 
model on scale free networks at temperature T = 0 as 
such studies have not been made so far to the best of our 
knowledge. Apart from the question whether the equilib¬ 
rium state is reached or not, we have also explored the na¬ 
ture of the state in case it does not. We are interested to 
see whether any active-absorbing phase transition occurs 
as the system parameters are varied. Both the random 
scale free network (RSF) and the Barabasi-Albert (BA) 
network have been considered for the study. Although 
many results are known for the WS network, we have ex¬ 
plored specifically the possibility of an active-absorbing 
phase transition in this network. The BA model is stud¬ 
ied to make direct comparison with the random scale free 
network, where results can be quite different [27]]. Also, 
comparison with respect to issues like freezing and ab¬ 
sorbing phase transition may be made for the RSF and 
the WS networks. 

In section II, we describe the network models and the 
dynamical evolution. The quantities which have been 
calculated are defined in Section III. The results are pre¬ 
sented in the next section and in the last section we sum¬ 
marise and discuss the studies made. 
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II. THE NETWORK MODELS AND DYNAMICS 

We have considered three different types of network: 
(a) Random scale free, (b)Barabasi-Albert and (3) Watts 
Strogatz (addition type) 3 network. We describe in brief 
how these networks are generated and the dynamical evo¬ 
lution process. In this section we also include a brief dis¬ 
cussion of how the numerical data had been analysed in 
previous studies. 

A. Random scale free (RSF) network 

In the random scale free network the degree distribu¬ 
tion follows a power law but otherwise the network is 
random. To generate random scale free network p7l - l32j . 
we assigned the degree of each node using the power law: 

V(k) ~ k~\ (1) 

where k is the degree of node and 7 is the characteris¬ 
tic degree exponent. The minimum value of k is 1 and 
the maximum cut-off value is y/N, where N number of 
nodes. This cut-off value ensures that there is no correla¬ 
tion [33| . Starting from the node with the maximum de¬ 
gree, links have been established with randomly selected 
distinct nodes. 

B. Barabasi-Albert (BA) network 

Barabasi-Albert network is a growing network where 
new nodes are joined to existing nodes with preferen¬ 
tial attachment. We start with the three fully connected 
nodes. Subsequently, a single node is added at a time 
to the network which is linked to one existing node. The 
probability that the new node is connected to the existing 
ith node with degree hi is given by [13, HH , 

n(fcj) = kj/y^kj. (2) 
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Degree distribution of the network is a power law with 
exponent 7 = 3; P{k ) ~ k~ 3 (27| . 

C. Watts Strogatz (addition type) model 

Addition type WS network is a one dimensional regu¬ 
lar chain with two nearest neighbour links as well as with 
some extra randomly connected long range links. Here 
the long range links have been added with probability 
q/N (total long range links ~ O(N)), where N is the 
number of nodes and q is a parameter, which denotes the 
number of extra long range links per node on an aver¬ 
age. So the average degree per node of this network is 
2 + q, which is a finite quantity as q/N —> 0 in the ther¬ 
modynamic limit. It is known Q that an order-disorder 
transition occurs at q = 1 . 


D. Dynamics of Ising model on networks 

The Hamiltonian of the Ising system in these networks 
can be expressed as 

H = — y ] JijSiSj, ( 3 ) 

i<j 

where S) = ±1 and Jij = 1 when sites i and j are con¬ 
nected and zero otherwise. Starting with the random 
configuration, single spin flip energy minimizing Glauber 
dynamics has been used to update the spin. In this dy¬ 
namics a randomly selected spin is flipped if the energy 
of the updated configuration is lowered. It is flipped with 
probability 1/2 if the energy remains unchanged on flip¬ 
ping. Fifty different network configurations have been 
considered and for each network hundred different initial 
configurations have been taken. The results are averaged 
over these configurations. We have considered system 
sizes N upto 1600. Periodic boundary condition has been 
used for the WS model which is embedded in real space. 

Previous studies have shown that for the thermally 
driven phase transition, not only mean field critical be¬ 
haviour exists for the Ising model on small world net¬ 
works but finite size scaling is also valid there IMH- The 
data collapse was obtained by rescaling the data and the 
scaling argument occurring in the scaling function was 
found to be of the familiar form eN 1 ^ (with u = 2) 
where e denotes the deviation from the critical point, v 
is argued to be equal to vd where v is the correlation 
length exponent and d the effective dimension of the sys¬ 
tem. With mean field critical exponent v = 0.5, d turns 
out to be equal to the upper critical dimension; d = 4. 

III. QUANTITIES CALCULATED 

We have estimated the following quantities in the 
present work. 

1. Magnetisation: m(t) has been calculated by taking 
the average of the absolute values of the magnetisation, 
m = (| ]Tb Si\)/N, as a function of time. Since evolution 
to both up spin dominated and down spin dominated 
configurations are possible, the absolute value is taken 
to compute the configuration average. Distribution of 
saturation value of magnetisation has also been estimated 
for random scale free (RSF) network. 

2. Residual energy: E r (t) = E(t) — E g where E g is the 
equilibrium energy of the ground state per spin and E(t) 
the energy per spin at time t. Residual energy measure¬ 
ment indicates the closeness to the equilibrium ground 
state. Since we are employing a zero temperature quench, 
the equilibrium ground state configuration corresponds 
to either all spins up or down. 

3- Pfii P (t) = Nfu p (t)/N, where Nfu p is the number 
of spin flips at time t, has been studied as a function of 
time. We count all the spin flips, i.e., if the same spin 
flips more than one time, all these occurrences are taken 
into account. 
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4. Freezing probability F is defined as the probability 
that the system does not reach the true ground state. 
It has been calculated as a function of the relevant pa¬ 
rameters. For T = 0, the known ground state is the 
ferromagnetic state with m(t —» oo) = 1. A frozen con¬ 
figuration will have an absolute value of magnetisation 
less than 1. 

5. Whether or not the dynamics continue, the system 
always reaches an iso-energy state. Time r to reach the 
iso-energy state has been calculated as a function of the 
relevant parameter for random scale free (RSF) network. 

The saturation values have been denoted using appro¬ 
priate subscript, e.g. m sa t for magnetisation. 


IV. RESULTS 


A. Random scale free network 
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FIG. 1. (Color online) RSF: Variation of freezing probability 
F( 7 ) as a function of 7 for different system size N. Inset 
shows the data collapse where F(j) has been plotted against 
A = ( 7 -7 l)N 1/0 . 

We first discuss the behaviour of the freezing probabil¬ 
ity F on the RSF. The variation of F with the network 
parameter 7 for different system sizes suggests that a 
freezing transition takes place here as the freezing prob¬ 
ability is unity above a value of 7 which we denote by 7 / 
(Fig. G}. Of course, if the initial state is fully ordered 
it is not counted as a frozen state. However, barring 
these two states (all up/ all down), none of the other ini¬ 
tial states reaches the actual equilibrium configuration. 
Hence the freezing probability is actually 1 — 2/2 N which 
becomes unity in the thermodynamic limit. Below 7 /, 
the freezing probability decreases with system size and 
shows a system size independent behaviour for 7 < 1 . 0 . 
The freezing probability is non zero for all values of 7 . 

As the freezing probability shows finite size dependence 
close to 7 / and since it is dimensionless, we argue that it 
should show a finite size scaling behavior in the following 
manner 


(4) 



m sat m sat 
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FIG. 2. (Color online) RSF: ((a), (b)) The distribution of 
saturation value of magnetisation for N = 576 for different 
values of 7 . (c) The variation of fluctuation of saturation 

value of magnetisation and (d) time to reach the energetically 
stable state as a function of 7 for different system size N. 
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FIG. 3. (Color online) RSF: Variation of magnetisation m{t) 
with time for different system size N. (a) 7 = 1.5 and (b) 
7 = 3.2. 


where gi is a scaling function. We estimate 7 / ~ 2.0 and 
the exponent v ~ 2.20 (inset of FigJT]). The estimates are 
obtained from the manifestly best collapse of the data 
after rescaling is done properly. As already discussed in 
section Hi D1 finite size scaling is valid for the Ising model 
on networks with v = vd where v is the mean field value 
of the correlation length exponent and d is equal to 4. 


F’(7,A r )=5i[(7-7 c / )^ 1/ 1 
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Assuming the same to hold good here, we get v k. 0.55 
from the fact that v ~ 2.2 which is fairly close to the 
mean field value 0.5 for Ising model [341] . 

We next show that an order-disorder phase transition 
also occurs in this network. The distribution of mag¬ 
netisation tt(m sat ) has been studied (Fig. 0r, b). It is 
unimodal in nature for the values of 7 far from the 7 /. 
The peaks occur at m sa t ~ 1 for 7 <C 7 / and m sa t ~ 0 
for 7 > 7 ^. The distribution is bimodal close to 7 / with 
the peak values at m sa t ~ 0 and m sa t ~ 1 . 

Fluctuations of the saturation value of the magneti¬ 
sation (Arri sa t 2 ) shows a peak which shifts toward 7 / 
and the the peak value increases as the system size is in¬ 
creased (Fig. [5}’). We have also estimated the time r to 
reach the energetically stable state. It shows a peak close 
to 7 £. Clearly r diverges near 7 ^, as the peak value in¬ 
creases with system size (Fig. 01). These behaviour sug¬ 
gest that there is an order-disorder transition also taking 
place as 7 is varied. In principle the order-disorder tran¬ 
sition may take place at a value 7 ™ ^ 7 / and we make 
further analysis to estimate 7 ™ more accurately. We also 
check whether a mean field-like behaviour is present for 
the order-disorder transition. 

The variation with time of the magnetisation m(t) has 
been shown for different system sizes for two different val¬ 
ues of the degree exponent 7 (Fig. 0i, b). The behaviour 
of the saturation values of the magnetisation m sa t for fi¬ 
nite sizes shows the typical characteristics of a continuous 
phase transition with 7 acting as the driving parameter 
(Fig. 0. 

Using finite size scaling method we indeed obtain a 
collapse of the data points for m sa t for different system 
sizes. The following scaling form for m sa t has been used: 

m sat = N-^g 2 [( 1 -^)N 1 /% (5) 

We obtain 7 ™ « 2.20, /3 ~ 0.47 and v ~ 2.20 (inset of 
Fig. 0. The values of the exponents /? and v are fairly 
close to the mean field values once again. 


and the spin flip probability which also attain a satura¬ 
tion value in time ( P sat and E sat ). Both the saturation 
values P S at and E sa t show nonmonotonic variation with 
7 (Figs. 0 0, with a peak which shifts as the system 
size is increased. P sa t in particular shows a very interest¬ 
ing behaviour with finite size. For 7 < 2.0, it decreases 
with system size indicating an absorbing phase which is 
actually the ordered phase as indicated by the behaviour 
of the magnetisation discussed above. However, there is 
a region between 7 « 2.0 and 7 ss 3.0, where P sa t in¬ 
creases with system size which indicates an active state. 
For 7 > 3.0, P sa t decreases with system size indicating an 
absorbing state once again. Hence we have two absorbing 
phases separated by an active phase and two transitions 
as the parameter 7 is varied. This is reminiscent of two 
distinct transitions observed in opinion dynamics models 


We found from the above studies that 7 / « 2.0 and 
7 ™ ~ 2.2 are close but not exactly the same. However, 
this could be due to finite size effects and these two might 
turn out to be identical, i.e., 7 ^ = 7 ™ = 2.0 in the ther¬ 
modynamic limit. This possibility is supported by the be¬ 
haviour of both r and E sat . While r diverges at 7 ~ 2.0, 
E sat decreases with system size below 7 ~ 2.0 (indicating 
an order phase) and increases or remains constant above 
this value (Figs. 0. 
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FIG. 4. (Color online) RSF: Variation of saturation value of 
magnetisation as a function of 7 for different system size N. 
Inset shows the data collapse where Y = m 3a tN^A has been 
plotted against X = (7 — 7 ™)V 1 / I/ . 

Next we discuss the behaviour of the residual energy 


FIG. 5. (Color online) RSF: Variation of saturation value of 
the fraction of spin flips as a function of 7 for different system 
size N. The vertical arrows separate different phases. 


B. Barabasi-Albert network 

The BA network has no intrinsic parameter. We calcu¬ 
late the relevant dynamic quantities for different system 
sizes. The magnetisation m(t) slowly increases for the 
first few time steps and gets saturated at long times for 
all the system sizes. Both E r (t) and Pfu p {t) decrease 
with time and then reach a saturation value. Saturation 
values of all these quantities as a function of system size 
have been plotted in Fig. 0i. The saturation value of 
magnetisation m sa t{N) decreases with system size and 
the variation shows a power law behaviour. We fitted 
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FIG. 6 . (Color online) RSF: Variation of saturation value of 
the residual energy as a function of 7 for different system size 


N. 


the variation with the form m sa t{N ) = a m N~ brn and the 
estimated exponents are a m ~ 1.467 and b m ss 0.150. 

We have also studied the variation of saturation value 
of the fraction of spin flips P sa t (N) and saturation value 
of residual energy E sat (N). Both of them show a slowly 
increasing behaviour with system size and these increas¬ 
ing behaviour also follow power law (Fig. |T}fi) . We fitted 
the variation of P sa t(N) with the form P sa t(N ) = afN bf 
and the exponents are a/ ~ 0.001 and bf ~ 0.346. For 
Esat(N), the fitted form is E sat (N) = aEN bE and the 
exponents are ob « 0.083 and ~ 0.293. The impor¬ 
tant point to note here is this seems to be an active phase 
as the saturation value of the fraction of spin flips shows 
increase with system size. In contrast, the same quanti¬ 
ties plotted for the RSF with 7 = 3.0 (Fig. |SJ) shows that 
it is an absorbing phase (as already noted in the previous 
subsection). We will discuss more about this observation 
in Section [V] The variation of E sat (N) with N in RSF is 
however, similar, with E sat (N ) = 0.0137V 0 279 . m sa t(N ) 
for RSF decays as 1.177V -0 408 which signifies a faster 
decay compared to the BA network. . 

In the BA model, the freezing probability F(N) in¬ 
creases with system size (Fig. [TJd). Freezing proba¬ 
bility —> 1 for N —> 00 . Thus the BA model is in 
an active disordered state where none of the configura¬ 
tion reaches the equilibrium ground state. The varia¬ 
tion of F(N) with system size fits well with the form 
F(N ) = 1 — e~ adN d and the calculated exponents are 
ad ~ 0.113 and bd ~ 0.548. The behaviour of the freez¬ 
ing probability as a function of N is also quite different 
for the BA and the RSF networks (with 7 = 3), in the 
latter we found a system size independent behaviour. 


C. WS network 

On the WS network, all the relevant quantities show 
a saturation behaviour as already noted previously on a 
slightly different version of the WS network [26|. 

Here in addition, we have studied the spin flip proba- 
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FIG. 7. (Color online) BA: (a) Variation of saturation value 
of magnetisation m sa t(N), fraction of spin flips P sa t{N) and 
residual energy E aa t(N) as a function system size, (b) Varia¬ 
tion of freezing probability F(N) as a function of system size 
N. 



FIG. 8 . (Color online) RSF (7 = 3.0): Variation of saturation 
value of magnetisation m sa t{N), fraction of spin flips P sa t(N) 
and residual energy E aa t(N) as a function system size. 


bilities. The saturation value of the probability of spin 
flips P sa t has been plotted against the parameter q for 
different system sizes and this quantity reveals interest¬ 
ing behaviour (Fig. Or). It is clearly seen that above 
q ~ 1, the probability decreases as a function of system 
size N while below q — l it is almost size independent. 
This indicates that there is an active absorbing phase 
transition taking place at this point. 

The value of the freezing probability F(q) either in¬ 
creases with N or remains constant which clearly shows 
that the entire phase is frozen for any q > 0. We find 
that for q < 1 the freezing probability reaches unity in 
the thermodynamic limit while it remains fairly constant 
beyond this value (Fig. [SJd). This is consistent with the 
active-absorbing phase transition stipulated to take place 
at q = 1; in the active state, one can never reach the equi¬ 
librium ground state configuration. 

F(q) shows a peak at q = 0.2 and the position of the 
peak is independent of system size which shows that the 
system is maximally disordered here. The peak values 
of freezing probability as a function of system size N 
is fitted with the form F pea k(N) = 1 — e~ aN and the 
estimated values of the exponents are a ~ 0.134 and b « 
0.540 (inset of Fig. [5 ]d). It may be noted that the same 
form is obeyed in the BA model. 
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FIG. 9. (Color online) WS: Variation of (a) saturation value 
of fraction of spin flips P sa t and (b) freezing probability F(q) 
as a function of q for different system sizes N (right panel). 
Inset of (b) shows the variation of peak value of F(q) with 
the system size. 
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V. DISCUSSIONS AND CONCLUSIONS 


We have studied zero temperature Glauber dynamics 
of the Ising model on three types of networks and com¬ 
pared the results. Frozen state is observed in all the three 
types of network models. For random scale free network, 
freezing probability is unity for 7 > 2 , i.e., the system 
never reaches the global equilibrium but for lower value 
of the parameter 7 , it decreases with system size and 
shows a system size independent behaviour for 7 < 1 . 
This behaviour of freezing probability suggests a freezing 
transition point at 7 / ~ 2.0. We also find an order disor¬ 
der transition point taking place very close to this point; 
in fact we believe that they occur at the same point and 
the difference is only a finite size effect. Also close to this 
point, the first active-absorbing (A-A) phase transition 
takes place; the disordered phase for 2 < 7 < 3 is active 
while for 7 < 2, one gets an absorbing phase. A second 
A-A transition takes place close to 7 ~ 3.0 and the sys¬ 
tem evolves to an absorbing disordered state beyond this 
value. In all probability these two A-A transitions take 
place at 7 = 2.0 and 3.0 in absence of finite size effects; 
these two points are significant as for 7 < 2 . 0 , the aver¬ 
age degree diverges while for 7 < 3.0, the degree variance 
diverges in the thermodynamic limit. 

One can compare the results of RSF network and 
BA network for the same characteristic degree exponent 
7 = 3.0. The residual energy shows an increase with sys¬ 
tem size N in both cases in a power law manner with 
an exponent which is fairly close. Also, the saturation 
value of magnetisation decreases with N in both net¬ 
works in a power law manner, corresponding exponents 
are however quite different. The saturation energy and 
magnetisation behaviour are consistent with the fact that 
the state is disordered for both RSF and BA networks. 
However freezing probability for BA model and RSF net¬ 
work show different behaviour with system size. For RSF, 


the freezing probability shows a system size independent 
behaviour while for BA model, it has a nonlinear depen¬ 
dence. As far as the spin flip probability is concerned 
RSF (at 7 = 3.0) again differs from the BA network. 
The RSF network and BA model are intrinsically differ¬ 
ent, BA is a growing network generated with a particular 
strategy. There is no loop in this network. In the case 
of RSF, the structure is completely different, there may 
be loops. Although in a numerical study, the value of 7 
may not be exactly 3 in either the RSF or BA networks 
due to finite size effects, we believe that this cannot be 
the reason for the results being qualitatively different. In 
fact previous studies on RSF and BA have shown that 
characteristic features may be quite different for the two 
networks HSHi. 

The results obtained for the WS model can also be 
compared to those found for the RSF network and BA 
network. In the WS network an A-A transition is ob¬ 
served at q ~ 1.0 which is also the order-disorder tran¬ 
sition point. Hence this is similar to the occurrence of 
an A-A and an order-disorder transition occurring simul¬ 
taneously in the RSF. However in the WS network, the 
entire disordered phase is active. For small average de¬ 
gree the system is in an active and for large degree in an 
absorbing phase. For the RSF network, two A-A tran¬ 
sitions exist where absorbing phase is observed for both 
large degree and small degree and an active state exists in 
between these two absorbing phases. In both RSF and 
WS another interesting feature is present, the ordered 
state shows a finite freezing probability with negligible 
system size dependence. In the WS network maximum 
value of the freezing probability F{q ) occurs at q ss 0.2 
and shows a behaviour similar to the freezing probability 
in BA network as a function of system size. 

To summarise, a systematic study of ordering dynam¬ 
ics of the Ising model on scale free networks has been 
made for the first time to the best of our knowledge. 
It is observed that the system freezes to a non equilib¬ 
rium steady state for all values of the relevant param¬ 
eters in the random scale free networks (RSF) and the 
Barabasi-Albert model (BA). The presence of two active¬ 
absorbing phase transitions in the RSF makes it differ¬ 
ent from the WS network where only one such transi¬ 
tion is observed. It is also concluded that in RSF one of 
the active-absorbing phase transition takes place at the 
order-disorder transition point which is similar to what 
is observed in the WS network. 
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